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Dissipative electromagnetic response and scattering geometry are potential sources for the ap¬ 
pearance of a negative Casimir entropy. We show that the dissipative contribution familiar from the 
plane-plane geometry appears also in the plane-sphere and the sphere-sphere geometries and adds to 
the negative Casimir entropy known to exist in these geometries even for perfectly reflecting objects. 

Taking the sphere-sphere geometry as an example, we carry out a scattering-channel analysis which 
allows to distinguish between the contributions of different polarizations. We demonstrate that 
dissipation and geometry share a common feature making possible negative values of the Casimir 
entropy. In both cases there exists a scattering channel whose contribution to the Casimir free 
energy vanishes in the high-temperature limit. While the mode-mixing channel is associated with 
the geometric origin, the transverse electric channel is associated with the dissipative origin of the 
negative Casimir entropy. By going beyond the Rayleigh limit, we find even for large distances that 
negative Casimir entropies can occur also for Drude-type metals provided the dissipation strength 
is sufficiently small. 

PACS numbers: 42.50.Lc, 05.70.Ce, 03.65.Yz 


I. INTRODUCTION 

The specific heat of a system coupled infinitely weakly 
to a thermal heat bath always has to be positive. This 
statement derives directly from the fact that the specific 
heat is proportional to the variance of the system energy. 
By appropriately integrating the specific heat, it follows 
that also the entropy has to be positive. Nevertheless, 
there are a number of situations in which the existence of 
a negative specific heat or a negative entropy is discussed. 

Such a situation immediately arises when the coupling 
between the system and the heat bath is no longer in¬ 
finitely weak dug. Then, the specific heat of the system 
can be defined as the difference between the specific heat 
of system plus heat bath on the one hand and the spe¬ 
cific heat of the heat bath alone on the other hand. While 
there is no reason why the difference of two positive num¬ 
bers has to be positive, one would typically expect that 
coupling additional degrees of freedom will increase the 
specific heat. The difference of specific heats just intro¬ 
duced is therefore expected to be positive. The same 
reasoning applies to the entropy. 

A notable exception is the damped free particle, for 
which under appropriate circumstances, the specific heat 
can become negative HH3HS]. An interesting feature of 
this system consists in the fact that its properties de¬ 
pend on the dimensionless ratio ksT/Ivy where T is the 
temperature and 7 the damping constant while ks and H 
are the Boltzmann constant and the Planck constant, re¬ 
spectively. As a consequence and in contrast to common 
expectations, decreasing the damping constant 7 will ren¬ 
der the free particle more classical and shift the transition 


* stefan.umrat h@physik.uni-augsburg.de 
t gert.ingold4Sphysik.uni-augsburg.de 


to the quantum behavior to lower temperatures. A nega¬ 
tive specific heat in the sense just discussed is also found 
in certain systems in condensed-matter physics and re¬ 
lated fields [HUlO]. 

Another important case where negative entropies give 
rise to considerable interest is the Casimir effect which we 
will consider in more detail in the present paper. Even 
though experimentally the Casimir force is much more 
relevant than the Casimir entropy, the former contains a 
significant zero-temperature contribution while the latter 
allows to focus on the effect of thermal photons. 

Suppose that we consider the Casimir interaction be¬ 
tween two objects. With each of the objects as well as 
with the combination of the two objects we can associate 
an entropy as the difference between the entropy of the 
electromagnetic field in the presence of the object(s) and 
the entropy of the free electromagnetic field. The Casimir 
entropy is an interaction entropy which then is obtained 
by subtracting the entropies of the two individual objects 
from the entropy of the two objects. While individual en¬ 
tropies are not finite, the Casimir entropy yields a finite 
value. 

One potential source of a negative Casimir entropy is 
the coupling of the electromagnetic field to the electrons 
in the scattering objects. These electrons undergo a dissi¬ 
pative motion and therefore the electromagnetic response 
of the objects is characterized by a damping constant 7 . 
The presence of dissipation inside the scatterers leads 
to a suppressed reflectivity of transverse electric (TE) 
modes at low frequencies while the scattering of trans¬ 
verse magnetic (TM) modes is not affected significantly. 
As a consequence, in the plane-plane geometry a reduc¬ 
tion of the Casimir force EH and Casimir entropy by a 
factor of two is found at high temperatures. The possi¬ 
bility of a negative Casimir entropy has been the subject 
of an extensive debate over more than a decade [T12H25] . 
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A thorough understanding of the role of dissipation in 
the Casimir effect is of considerable interest for the inter¬ 
pretation of experiments. Presently, there still exists dis¬ 
agreement about whether the dissipative low-frequency 
part of the electromagnetic response is relevant or not. 
While some experiments are interpreted in terms of the 
dissipation-less plasma model ESI [32, other experiments 
are interpreted in terms of the Drude model which in¬ 
cludes dissipation [321 [TTTTj . At this point, it is useful to 
recall that a nonvanishing value of 7 is required for the 
Lifshitz formula to be applicable [33] , 

Most Casimir force measurements are carried out in 
the plane-sphere geometry with a surface-to-surface dis¬ 
tance between 0.1 and a few micrometers. Depending on 
the experiment, the radius of the sphere or the radius 
of curvature of a spherical segment ranges from 20 fj,m 
[33 SS] to about 20 cm E2- The latter case is rather close 
to a configuration of two parallel plates, where TE and 
TM modes do not mix. However, with decreasing sphere 
radius, polarization mixing becomes increasingly relevant 
and it is important to understand how the occurrence of 
negative Casimir entropies depends on dissipation in the 
plane-sphere and sphere-sphere geometry. 

The extension of the discussion of the negative Casimir 
entropy due to dissipation beyond the plane-plane geom¬ 
etry is not the only purpose of this paper. In fact, it is 
known since some time, that negative Casimir entropies 
can arise for the plane-sphere [221 [23] [25] and the sphere- 
sphere m geometry even for perfect conductors. An 
important question therefore concerns the relative status 
of the negative Casimir entropies due to dissipation on 
the one hand and due to geometry on the other hand. 

Recently, it was shown that the negative Casimir en¬ 
tropy in the perfect metal plane-sphere and sphere-sphere 
geometries can be traced back to the mixing between the 
TE and TM modes of the electromagnetic field [231[29] ■ 
This identification became possible by realizing that neg¬ 
ative Casimir entropies are most pronounced at large 
distances between the scattering objects. Then, within 
a scattering formalism, it is sufficient to consider one 
round-trip of electromagnetic waves between the scat¬ 
tered and a decomposition into the various scattering 
channels becomes possible. 

In Refs. [25] [25] it was argued that in the large-distance 
limit two Drude-type metal objects will not allow for a 
negative Casimir entropy because the reflection of the TE 
modes is strongly suppressed and only a positive contri¬ 
bution of TM modes is left. This regime can be attained 
for any given damping strength by making the distance 
between the scattering objects sufficiently large. On the 
other hand, for any given distance of the objects, one 
can make the damping strength so small that a negative 
Casimir entropy ensues. 

In order to disentangle the contributions to a negative 
Casimir entropy arising from dissipation and geometry, 
we will go beyond the electric dipole approximation, i.e. 
the Rayleigh limit [35] employed in Ref. [25]. While we 
will assume the distance d between the scattering objects 


to be large compared to the radii of the spherical objects 
involved, we will allow for damping strengths so small 
that a dissipative contribution to the negative Casimir 
entropy can occur in addition to the geometric contribu¬ 
tion. 

In Sec. [Tl] we start by sketching the formalism re¬ 
quired for the polarization channel analysis introduced 
in Ref. [25]. By considering the Casimir entropy for 
the plane-plane, sphere-plane, and sphere-sphere geome¬ 
tries, we then illustrate how negative Casimir entropies 
of dissipative and geometric origin manifest themselves 
in the temperature dependence of the Casimir entropy. 
In Sec. Ill the contributions from the various scattering 
channels are analyzed within the dipole approximation 
but beyond the Rayleigh limit. The approximation of a 
single scattering round-trip between the two scattering 
objects allows us to disentangle dissipation and geome¬ 
try as sources of a negative Casimir entropy and at the 
same time to identify common features. In Sec. m we go 
beyond the single round-trip approximation and demon¬ 
strate that repeated scattering round-trips and, for per¬ 
fect conductors, higher multipoles tend to suppress neg¬ 
ative Casimir entropies. In Sec. [V] we present our conclu¬ 
sions. An appendix collects some formulae needed for the 
evaluation of the Casimir entropy within the scattering 
formalism. 


II. DISSIPATIVE AND GEOMETRICAL 
INFLUENCE ON CASIMIR THERMODYNAMICS 


The Casimir entropy is obtained from the Casimir free 
energy T by means of the usual thermodynamic relation 


$ = - 


a f _ 

dT ' 


(1) 


Within the multiple-scattering theory [391442] . the 
Casimir free energy can be expressed as 

Jr=k ^~ E Tr (in [l — A4 (|£„ |)]) (2) 


with the Matsubara frequencies £„ = 2nnkBT/h. The 
matrix A4 describes round-trips of electromagnetic waves 
between all involved scatterers at imaginary frequencies 

i Sn ■ 

As we will focus on the sphere-sphere and plane-sphere 
geometries, the trace in ( 2 ]) can be taken in a spherical 
wave basis with the angular part given by spherical har¬ 
monics of degree l and order m. In view of the axial sym¬ 
metry of the geometries considered here, m is conserved 
during the scattering process. The round-trip operator 
in a corresponding subspace 


M (m) = 


(3) 


can be expressed as product of translation operators 77y 
from the reference frame of object j to that of object 
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i and reflection operators IZi associated with object i. 
For reference, explicit expressions for the sphere-sphere 
geometry involving Mie coefficients and spherical wave 
translation formulas are given in the appendix. 

For large distances between the scattering objects, all 
matrix elements of the round-trip operator M will be 
much smaller than one. The Casimir free energy § can 
then be linearized, 

U rp oo 00 

£ £ ™ (M) (I^I)> (4) 

n =—00 m =—00 

thereby retaining only contributions from single round- 
trips. Within this approximation it is possible to decom¬ 
pose the Casimir free energy T and the Casimir entropy S 
into contributions from the different scattering channels 
and thereby to gain physical insight j2U G2| ■ Of particu¬ 
lar interest will be the channel where the TE polarization 
is conserved during the round-trip and the polarization¬ 
mixing channel where each of the two objects scatters a 
different polarization. 

Furthermore, for sufficiently large separations d/R > 
20, the dipole approximation applies to the Casimir en¬ 
tropy and we can restrict ourselves to £= 1 and \m\ = 0,1 
(see Figs. [5] and [7] below). We will make use of the single 
round-trip dipole approximation in the present section as 
well as in the following Sec. m Its range of validity will 
be discussed in Sec. ITVl 

In order to obtain a first idea of the interplay between 
the dissipative and geometric contributions to the neg¬ 
ative Casimir entropy in the plane-sphere and sphere- 
sphere geometries, we compare with the plane-plane ge¬ 
ometry where the negative Casimir entropy is solely 
of dissipative origin. Apart from the special case of 
perfectly conducting metals (P), we will consider non¬ 
magnetic Drude-type (D) metals with the frequency- 
dependent permittivity 


e(w) = 1 - 


w(w + ij) 


(5) 


Here, 7 and ujp are the damping constant and the plasma 
frequency, respectively, corresponding to a dc conductiv¬ 
ity (To = Wp/ 7 . While up determines the frequency scale 
beyond which the metal becomes transparent, a finite 
value of (To results in a vanishing reflectivity of the TE 
modes at zero frequency. 

Figure [I] shows the temperature dependence of the 
Casimir entropy for the three geometries indicated by 
the insets. The entropy is scaled by its high-temperature 
value iSy T for perfectly conducting metals. A dimension¬ 
less temperature is defined by means of the distance d be¬ 
tween the reference frames of the two scattering objects. 
The three solid lines correspond to Drude-type metals 
with jd/c = 10~ 2 ,10 2 , and 10 4 increasing from bottom 
to top and ujpd/2nc = 400. Here, c is the speed of light. 
The lower (blue) curve and the upper (red) curve thus 
correspond to good and bad conductors, respectively. As 
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FIG. 1. The Casimir entropy S scaled by its high-temperature 
value SK t f° r perfect conductors is shown as a function of 
temperature for the (a) plane-plane geometry, (b) sphere- 
sphere geometry, and (c) plane-sphere geometry. While the 
dashed lines represent the entropy for perfectly conducting 
objects, the solid lines correspond to Drude-type conductors 
with the damping constant 'yd/c = 10 ~ 2 , 10 2 , and 10 4 increas¬ 
ing from bottom to top and the plasma frequency given by 
ujpd/2-KC = 400. The sphere radii are chosen as R = d/20. 
The arrows indicate the effect of the missing zero-frequency 
contribution to the Casimir entropy for Drude-type metals 
with respect to perfect conductors. 


a reference, the case of a perfectly conducting metal is 
depicted as dashed line. 

We start by recalling in Fig. o> the main features of 
the Casimir entropy for Drude-type metals in the plane- 
plane geometry. It is well known that a finite dc conduc¬ 
tivity suppresses the zero-frequency contribution of the 
TE modes, resulting in a reduction of the Casimir 

entropy at high temperatures by a factor of two with re¬ 
spect to the case of an infinite dc conductivity m The 
arrows in Fig. [l] visualize the effect of this missing zero- 
frequency contribution to the Casimir entropy. 

Since for perfectly conducting metals the Casimir en- 
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tropy decreases with decreasing temperature, the miss¬ 
ing zero-frequency contribution of the TE modes exceeds 
their contribution at finite temperatures. As a conse¬ 
quence, at sufficiently low temperatures, the Casimir en¬ 
tropy for Drude-type metals becomes negative, provided 
the clc conductivity is sufficiently large. For the smallest 
value of the damping constant 7 in Fig. the entropy 
curve over a large temperature range ressembles closely 
that for a perfect conductor except for the entropy shift 

A^nf • 

The Casimir entropy approaches its vanishing zero- 
temperature value only when the temperature falls be¬ 
low a temperature of the order of fty/fes [2U]. The mid¬ 
dle (black) curve therefore comes close to zero at much 
higher temperatures than the lower (blue) curve. This 
behavior is reminiscent of the damped free particle dis¬ 
cussed in the introduction where the temperature scale 
is also set by the damping constant 7 . 

A reduction of the Casimir entropy at high temper¬ 
atures as well as negative values at lower temperatures 
with a crossover to vanishing Casimir entropy on a tem¬ 
perature scale given by 7 are also visible for the sphere- 
sphere geometry in Fig. [Tjr and the plane-sphere geome¬ 
try in Fig. [^. Since the two geometries share the same 
qualitative features, we restrict the latter figure to neg¬ 
ative values of the Casimir entropy. In both figures, the 
sphere radii are chosen as R = d/20 to ensure the validity 
of the single round-trip dipole approximation. 

There are quantitative differences in the temperature 
dependence of the Casimir entropy for the three geome¬ 
tries. While the temperature scale on which the Casimir 
entropy approaches zero at low temperatures clearly in¬ 
creases with increasing damping constant 7 , its value 
depends on the scattering geometry. In addition, the 
reduction of the Casimir entropy at high temperatures 
amounts only to 1/3 for the plane-sphere geometry and 
1/5 for the sphere-sphere geometry instead of 1/2 for the 
plane-plane geometry. This reduction is a consequence of 
the fact that for a perfectly conducting sphere the reflec¬ 
tivity of the TE modes is only one half of the reflectivity 
of the TM modes. 

Besides the features already present in the plane-plane 
geometry, an additional dip in the Casimir entropy ap¬ 
pears for the sphere-sphere and plane-sphere geometries 
at temperatures where 27 t ksTd/hc is of order one. This 
feature is of geometric origin and survives even in the ab¬ 
sence of dissipation as can be seen from the dashed lines 
in Figs. [T] > and c. Its position depends on the geometry 
and it is more pronounced in the sphere-sphere geometry. 

While for a damping constant as large as 7 = 10 4 c/d, 
i.e. for the upper (red) solid curves, the Casimir entropy 
still becomes negative in the plane-sphere geometry, it re¬ 
mains positive for all temperatures in the sphere-sphere 
geometry. Nevertheless, even there a remnant of the dip 
is visible. However, for sufficiently large damping con¬ 
stant 7 , the contribution to the Casimir entropy involv¬ 
ing TE modes will disappear [281 ^5] , For such strong 
damping, the electrons will suffer many collisions during 


an oscillation period of the electromagnetic field for all 
relevant frequencies. The total irradiating power in the 
TE modes will then be dissipated, thereby suppressing 
the reflection of these modes. 


III. CHANNEL ANALYSIS BEYOND THE 
RAYLEIGH LIMIT 

The temperature dependence of the Casimir entropy 
displayed in Figs. [T]) and c suggests that the negative 
Casimir entropy can be separated into contributions of 
dissipative and of geometric origin. This is particularly 
clear for the sphere-sphere geometry for which we will 
demonstrate in the second half of this section that the 
expected separation is indeed possible in terms of differ¬ 
ent scattering channels. Even though we will be able to 
disentangle geometric and dissipative parts, Figs. |T] > and 
c also show that the geometric contribution depends on 
the damping constant appearing in the Drude-type per¬ 
mittivity ([ 5 ]). In order to analyze this dependence, we 
will need to take a closer look at the scattering proper¬ 
ties of Drude-type metal spheres. 

For the sphere-sphere geometry, the translation oper¬ 
ators T appearing in |3]) restrict the relevant wave num¬ 
bers to small values with kd < 1. This cut-off is cau sed b y 
the modified spherical Bessel functions present in (A.4|. 


For large distances d R, one may therefore apply the 
Rayleigh limit, where only dipole scattering is relevant. 
For perfect metal spheres, it will be sufficient to consider 
only the leading term in the expansion of the l = 1 Mie 
coefficients for TM modes 

a* = -^(kR) 3 + ^(kR) 5 + 0(k 6 ) , ( 6 ) 

as well as for TE modes 


= _( kR) 3 + ~{kR) 5 + O ( k 6 ) . 


(7) 


Here, TM and TE scattering amplitudes are of the same 
order with a relative factor of two alluded to in the pre¬ 
vious section. 

In contrast, for Drude-type metal spheres, the Mie co¬ 
efficients for TM modes with i=l 

af = -^(kR) 3 + ^(kR) A + 0{k 5 ) ( 8 ) 

and for TE modes 


1 

’ 2 

fa 0 R > 

2 

45 

21 

K~j 

+ 7 j 


{kR) 5 + 0{k 6 ), 


(9) 


differ in the power of kR in the leading term. The scatter¬ 
ing of TE modes thus becomes negligible in the Rayleigh 
limit. However, the prefactor of the leading term in 
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indicates the appearance of a new scale. As the compari¬ 
son of ^ and 0 shows, the scattering of TE modes can 
only be neglected if 


Ran d 

-- < — 

30c R 


( 10 ) 


The applicability of the Rayleigh limit thus also depends 
on the dc conductivity (To. For a given dc conductivity, it 
is always possible to reach the large-distance limit where 
the scattering of TE modes can be neglected. On the 
other hand, for a given distance, the Rayleigh limit ceases 
to hold if the dc conductivity becomes too large. 

A physical interpretation of the condition (10) can be 
given in terms of the quasi-static skin depth 



/ 2c 2 \ 1/2 

*M= — • (ii) 

\a 0 to J 

Scattering of TE modes can thus only be neglected pro¬ 
vided the sphere is sufficiently small compared to the skin 
depth, 


Vl5S, 


( 12 ) 


in the relevant frequency range. Then, no eddy currents 
can be induced and the sphere becomes transparent to 
the TE modes. 

The breakdown of the validity of the Rayleigh regime 
for a well conducting metal sphere with a^R/c = 47T x 
10 4 is illustrated in Fig. [ 2 ] The ratio of Mie coefficients 
—/af for dipole scattering approaches at small wave 
numbers the behavior predicted by the leading terms in 
| 8 ]) and ©• At larger wave numbers, the dipole regime 
is entered where the ratio of the Mie coefficients af and 
bf approaches the value for perfectly conducting spheres 
given by 0 and 0 - Increasing the wave numbers even 
further, the multipole regime is reached where first the 
scattering of TM modes with £ = 2 described by the Mie 
coefficient becomes relevant. 

An increase in the dc conductivity will reduce the 
Rayleigh regime and render the dipole regime valid at 
even smaller frequencies. To ensure the validity of the 
Rayleigh limit, one thus has to choose either bad con¬ 
ductors or go to large distances satisfying (10 1 . In the 


large-distance limit considered in Refs. [28] and [29], it 
was assumed that (10) holds and that the TE modes are 


not scattered by Drude-type metal spheres. Under this 
condition, it was possible to derive analytical results for 
the Casimir entropy. Here, we will allow for an arbitrary 
dc conductivity and thus have to resort to a numerical 
evaluation of the Casimir entropy employing the full Mie 
coefficients given in the appendix. 

In the large-distance limit d > R, where the linearized 
single round-trip expression Q for the Casimir free en¬ 
ergy applies, we can decompose the Casimir entropy into 
contributions arising from the scattering channels char¬ 
acterized by the polarization on the two spheres. For 
two channels the polarization on the spheres is conserved 


FIG. 2. The ratios of Mie coefficients — 6°/a? and a®/af 
as a function of the wave number k are shown for a Drude- 
type metal sphere with radius R. The plasma frequency 
and the damping constant are given by ujpR/2nc = 20 and 
'y/wp = 10 -4 . Three different regimes can be distinguished. 
In the Rayleigh regime, TE scattering is negligible, while in 
the dipole regime, TE and TM scattering are of comparable 
strength. In the multipole regime, the scattering of higher 
multipole waves becomes relevant. The dashed lines indicate 
the approximations valid in the first two regimes with refer¬ 
ences to the corresponding equations. 


along the round-trip: TM^TM and TE^TE. In the fol¬ 
lowing, we will refer to them as TM and TE channels, 
respectively. In addition, two channels correspond to po¬ 
larization mixing: TM^=»TE and TE^iTM. Because of 
the symmetry of our geometry with identical spheres, 
the two polarization-mixing channels provide equal con¬ 
tributions and we will only consider their sum. 

Fig. [3] displays the temperature dependence of the 
three different contributions to the Casimir entropy 
scaled by the high-temperature Casimir entropy <S|( T 
for perfectly conducting spheres. The parameters equal 
those used in Fig. [T| but no curves for the highest dc 
conductivity are shown. 

In Fig. [3jr, the contribution of the TM channel is 
depicted. We only show the curve corresponding to 
7 d/c = 10 2 but omit the curves for 'yd/c = 10 4 and 
for perfect conductors. According to ^ and the 
reflection coefficient for TM modes depends only very 
weakly on the dc conductivity and, in fact, the other 
curves would lie within not more than one line width. As 
discussed before, at high temperatures the TM channel 
contributes 4/5 of the total Casimir entropy. The verti¬ 
cal axes in Fig. [3] are drawn to scale so that a comparison 
between the three contributions is facilitated. While the 
TM channel largely dominates at high temperatures, this 
is no longer the case at low temperatures, thus opening 
the possibility for negative Casimir entropies. 

The temperature dependence of the TE channel shown 
in Fig. 0 > displays the same features as the Casimir en¬ 
tropy in the plane-plane geometry depicted in Fig. 0 
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FIG. 3. The contributions from (a) the TM channel, (b) the 
TE channel, and (c) the polarization-mixing channels to the 
Casimir entropy are shown as function of temperature for a 
sphere-sphere geometry with d/R = 20. All contributions are 
scaled by the total high-temperature Casimir entropy Sh T 
for perfectly conducting spheres. The parameters are cho¬ 
sen as in Fig. |Tj i.e. the dashed lines correspond to perfect 
metal spheres while the solid black and red lines correspond 
to Drude-type metal spheres with 7 d/c = 10 2 and 10 4 , respec¬ 
tively, and ujpd/2nc = 400. For the sake of clarity, panel (a) 
displays only one line because all other lines deviate by not 
more than one line width. To facilitate the comparison be¬ 
tween the contributions of the three types of scattering chan¬ 
nels, the scales on the vertical axes are chosen to be equal. 


once the contribution arising from the TM modes is sub¬ 
tracted. We recognize the shift between the Casimir en¬ 
tropies for perfectly conducting spheres and Drude-type 
metal spheres at high temperatures due to the vanish¬ 
ing reflectivity for the TE modes at zero frequency. For 
not too small dc conductivity, i.e. the solid black curve, 
the contribution to the Casimir entropy becomes negative 
and goes to zero on a temperature scale proportional to 
7 . For bad conductors, i.e. the red solid curve, the contri¬ 
bution of the TE channel almost vanishes, in agreement 
with our discussion of the Rayleigh limit. The scattering 
channel analysis thus confirms that the TE modes are 
not only responsible for the dissipative contribution to 
the negative Casimir entropy in the plane-plane geome¬ 
try, but also in the sphere-sphere geometry. It is worth 
mentioning that the same conclusion can be shown to 
hold true for the plane-sphere geometry. 

Polarization-mixing channels do not exist in the plane- 


plane geometry but contribute to the Casimir entropy 
in the plane-sphere and sphere-sphere geometries. The 
temperature dependence for the latter case is shown in 
Fig#- From Refs. J25] and [29] it is already known that 
the polarization-mixing channels are associated with the 
geometric origin of a negative Casimir entropy. Here, 
we can now study the dependence of this feature as a 
function of the Drude damping constant 7 . 

As the curves in Fig. Hf 5 show, the negative Casimir 
entropy tends to be increasingly suppressed with increas¬ 
ing damping constant 7 . In constrast to the TE channel 
in Fig. S>, a continuous transition exists from perfectly 
conducting spheres, where the negative Casimir entropy 
is most pronounced, to bad Drude-type metal spheres, 
where the mode mixing is significantly suppressed. Re¬ 
ducing the dc conductivity even further would ultimately 
result in a vanishing geometric contribution to a negative 
Casimir entropy [SSJ US]. 

We conclude that the geometric and dissipative contri¬ 
butions to a negative Casimir entropy can clearly be dis¬ 
tinguished by means of a scattering-channel analysis pro¬ 
vided the single round-trip condition is met, i.e. d > R. 
Nevertheless, the two types of channels share a common 
feature as can be seen from Fig. [4] where the contribu¬ 
tions of the different scattering channels to the Casimir 
free energy is shown as a function of temperature. 

The contributions of both, TE and polarization-mixing 
channels vanish in the high-temperature limit. This be¬ 
havior is in clear contrast to the usual behavior found for 
the TM channel which yields a negative contribution de¬ 
creasing linearly with temperature at high temperatures. 
In view of # and the fact that the contributions to the 
Casimir free energy at zero temperature are negative, the 
temperature dependence of the first two channels nec¬ 
essarily implies a negative contribution to the Casimir 
entropy in some temperature range. 

This common feature should not obscure their differ¬ 
ent physical origins. For the TE channel, the contribu¬ 
tion to the Casimir entropy at high temperatures van¬ 
ishes because of the vanishing zero-frequency reflectivity. 
In contrast, for the polarization-mixing channel, no con¬ 
tribution to the Casimir entropy can occur at zero fre¬ 
quency where the fields corresponding to the two modes 
are either purely electric or purely magnetic in any of 
the two coordinate systems centered at each sphere. The 
difference in the physical mechanism manifests itself in 
Figs. and c, where we observe that the transition from 
a Drude-type metal to a perfect conductor occurs contin¬ 
uously for the polarization-mixing channel while it does 
not for the TE channel. 


IV. SPHERE-SPHERE CASIMIR ENTROPY 
BEYOND THE DIPOLE APPROXIMATION 

So far, our discussion of the Casimir entropy has relied 
on the large-distance assumption d R. Only then it 
suffices to consider single round-trips and the contribu- 












7 



FIG. 4. The contributions of the different scattering chan¬ 
nels to the scaled free energy T = (d/RfT are shown as 
a function of temperature. The distance between the two 
Drude-type metal spheres is given by d/R = 20 and the ma¬ 
terial parameters are yd/c = 10 and ojpd/2nc = 400. The 
inset displays the low-temperature behavior of the TE chan¬ 
nel and demonstrates that the corresponding contribution to 
the Casimir entropy vanishes in the zero-temperature limit. 

tions of dissipative and geometric origin can be added. 
When the distance between the two spheres is decreased, 
it becomes necessary to go beyond the dipole approxima¬ 
tion and to allow for scattering processes involving more 
than a single round-trip. 

In this section, we first discuss the role of multiple 
round-trip contributions for perfectly conducting spheres 
and thus concentrate on the geometrically induced nega¬ 
tive Casimir entropy. In a second step, we add dissipation 
by considering Drude-type metal spheres. 

Fig. [5] shows the Casimir entropy as a function of 
the distance d between two perfectly conducting spheres. 
The temperature is kept fixed at 2iTkBTR/hc = 1. The 
black points correspond to the full calculation based on 
([2]) including multipoles up to £ ma x = 60 while the grey 
points correspond to the single round-trip approximation 
©• The white points are obtained by additionally ap¬ 
plying the dipole approximation. 

The inset of Fig. [5] shows how the single round-trip 
dipole approximation approaches the exact result for the 
Casimir entropy as the distance between the two spheres 
is increased. At d/R = 20, the distance chosen in the 
previous sections, the relative error reaches about one 
percent. 

At the temperature chosen here, the Casimir entropy is 
positive at large distances between the two spheres. De¬ 
creasing the distance, the entropy will eventually become 
negative and go through a minimum before rising again 
to positive values m- Such a distance dependence was 
also found for the plane-sphere geometry [22J [23]. Inter¬ 
estingly, even the single round-trip dipole approximation 
(white points) is capable of qualitatively describing this 
distance dependence of the Casimir entropy. 

Including higher multipoles and multiple round-trips 



FIG. 5. The Casimir entropy S is presented as a function of 
the distance d between two perfectly conducting spheres of 
radius R at 2nksTR/hc = 1 for different levels of approxima¬ 
tion. Results obtained within the single round-trip approx¬ 
imation are displayed as white and grey points. While the 
first ones correspond to the dipole approximation, the latter 
include higher multipoles up to £ max = 60. The numerically 
exact result depicted by black points is obtained on the basis 
of 0 with £ max = 60. The inset shows the relative difference 
of the dipole approximation S dlp and the exact result S ex . 

yields significant positive contributions to the Casimir en¬ 
tropy in particular at relatively small distances between 
the two spheres. However, already the single round-trip 
approximation (grey points) provides a good description 
almost down to distances where the minimum of the 
Casimir entropy is reached. Only at even smaller dis¬ 
tances multiple round-trip contributions become relevant 
(black points). 

Fig. 0 clearly demonstrates that for not too small 
distances d the main correction to the single round- 
trip dipole approximation consists in the contribution of 
higher multipoles. We therefore analyze in Fig. [6] how 
the dipole-quadrupole and quadrupole-quadrupole chan¬ 
nels contribute to the temperature dependence of the 
Casimir entropy. The distance between the two perfectly 
conducting spheres is chosen as d = 2.75 R, i.e. close to 
the minimum of the Casimir entropy in Fig. [5] 

The solid line depicts the Casimir entropy including 
multipoles up to £ max = 15 while the contributions 
due to dipole-dipole, dipole-quadrupole, and quadrupole- 
quadrupole scattering are displayed by the dashed, dash- 
dotted, and dotted lines, respectively, and comprise all 
polarization channels. In agreement with Ref. [28] we 
find that the negative Casimir entropy is caused exclu¬ 
sively by dipole scattering. In fact, as can be seen from 
Fig- i for dipole-quadrupole and quadrupole-quadrupole 
scattering, higher multipoles lead to positive Casimir en¬ 
tropy contributions. Still, an analysis of the polariza¬ 
tion channels involved in the multipole contributions de¬ 
picted in Fig. [6] reveals that they have the same signs as 
the corresponding dipole-dipole contributions. However, 
for higher multipoles the negative contribution from the 
polarization-mixing channel is shifted to larger temper- 
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FIG. 6. The scaled Casimir entropy S = ( d/R) 6 S as a func¬ 
tion of temperature is shown as solid line accounting for multi¬ 
poles up to tax = 15. The contributions arising from dipole- 
dipole, dipole-quadrupole, and quadrupole-quadrupole scat¬ 
tering are depicted as dashed, dash-dotted, and dotted line, 
respectively. The distance between the perfectly conducting 
spheres is d = 2.75 R. 


atures where the polarization-conserving channels domi¬ 
nate, yielding an overall result with positive sign. 

We now turn to the dissipative case and consider 
Drude-type metal spheres. The distance dependence 
of the Casimir entropy is depicted in Fig. [7] for 
2nk B TR/hc = 1. The material parameters correspond 
to those of the best conductor considered in Sec. El and 
represented by the blue line in Fig. jTj Only results 
within the single round-trip approximation are displayed 
in Fig. [TJ with white and grey points representing, respec¬ 
tively, the dipole approximation and the Casimir entropy 
with higher multipoles included. For the data presented 
here, the numerically exact result accounting for multiple 
round-trips practically coincides with the grey points. 

In contrast to perfectly conducting spheres, the contri¬ 
butions arising from higher multipoles are now negative 
and lead to a Casimir entropy increasing monotonically 
with distance. While the geometrically induced negative 
Casimir entropy is only due to dipole scattering, this is 
not the case for the negative Casimir entropy caused by 
dissipation. Here, also higher multipole scattering con¬ 
tributes with negative values. This important difference 
is to be expected because dissipation can give rise to a 
negative Casimir entropy even in the plane-plane geom¬ 
etry. 

As the inset in Fig. [7] shows, the deviation of the sin¬ 
gle round-trip dipole approximation from the numerically 
exact result is comparable for, both, perfect conductor 
spheres and Drude-type metal spheres. However, because 
of the missing high-temperature contribution of the TE 
modes, the absolute value of the Casimir entropy in the 
dissipative case exceeds the one for perfect conductors by 
one order of magnitude for the temperatures considered 
here. As a consequence, multiple round-trip corrections 
are much less important in Fig.[7]than they were in Fig. [5] 


FIG. 7. The Casimir entropy S is presented as a function 
of the distance d between Drude-type metal spheres of ra¬ 
dius R at 2-nksTR/hc = 1 within the single round-trip ap¬ 
proximation. The material parameters correspond to those 
of the best conductor represented by a blue curve in Fig. [I] 
i.e. ojpR/2nc = 20 and ujp/2n-y = 4 x 10 4 . The white points 
correspond to the dipole approximation while the grey points 
include higher multipoles up to f m ax = 60. Numerically exact 
results accounting for multiple round-trips cannot be distin¬ 
guished from the grey points in this figure and are therefore 
not shown. The inset displays the relative difference of the 
dipole approximation <S dlp and the exact result S ex . 


However, in both cases, the approximations central to 
our analysis - single round-trip and dipole scattering - 
are well justified for the large distances d/R = 20 chosen 
in the previous sections. 


V. CONCLUSIONS 

Dissipation and scattering geometry represent possi¬ 
ble sources of negative Casimir entropies. While the first 
mechanism has been studied for quite some time in the 
plane-plane geometry, the second mechanism was identi¬ 
fied for plane-sphere and sphere-sphere geometries with 
perfect conductors. For Drude-type metals of sufficiently 
large clc conductivity, we have shown the coexistence of 
both mechanisms by going beyond the Rayleigh limit. 

In the large-distance limit, where only single round- 
trips are relevant, we employed a scattering-channel anal¬ 
ysis to disentangle dissipative and geometric origins of 
the negative Casimir entropy. Concentrating on the 
sphere-sphere geometry, we found that dissipation acts 
via the TE channel as in the plane-plane geometry. In 
contrast, the geometric contribution can be traced back 
to the polarization-mixing channel as in the perfectly 
conducting case. Although being of different physical 
origin, both mechanisms are associated with a vanishing 
free energy in the high-temperature limit. 

Focussing on the geometrical origin of negative Casimir 
entropies, we have found for perfectly conducting scat- 
terers that while the dipole-dipole channel can lead to 
















9 


a negative contribution to the Casimir entropy, higher 
multipoles yield a positive contribution. With decreasing 
distance between the scatterers, the higher multipoles be¬ 
come dominant and the total Casimir entropy turns posi¬ 
tive, thereby making connection with the plane-plane ge¬ 
ometry for perfect conductors. In contrast, for the dissi¬ 
pative mechanism leading to negative Casimir entropies, 
the negative contribution is not restricted to dipole scat¬ 
tering. 
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Appendix: Reflection and translation operators 


For the reflection at a sphere, the scattering opera¬ 
tor 1Z occurring in |3| is naturally represented in the 
spherical multipole basis characterized by t and to. The 
corresponding matrix elements are given by the Mie co¬ 
efficients (58) . At imaginary frequencies and for non¬ 
magnetic spheres, p = 1, with refractive index n = y/e, 
they can be expressed as 


. . e+1 n nip^nkRjip'^kR) — ipi(kR)ij/ t (nkR) 

1 2 nipt(nkR)£,' e (kR) — kR)ip , i (nkR) 

(A.l) 

for TM polarization 


i 1 ^+i 7r ipe(nkR)^' e (kR) — nipe{kR)'ilj' e (nkR) 
e 2 il>e(nkR)£,' i (kR) — n^£(kR)ijj , e (nkR) 

(A.2) 

for TE polarization. The Riccati-Bessel functions are 
defined in terms of modified spherical Bessel functions of 
the first kind and second kind [55] as: 


ipe{p) = pie(p), &(p) = pk e (p). (A.3) 


The spherical vector wave basis refers to a coordinate 
system with origin in the center of one of the spheres. 
A basis change from one center to another one can be 
done by means of translation formulae [55H571 which can 
also be expressed in terms of imaginary wave vectors |42| . 


The matrix elements of the translation operator take a 
relatively simple form if the translation occurs along the 
z axis. Then, to is conserved while t and the polarization 
P can change. For imaginary wave numbers fc, the matrix 
elements can be cast into the form 


T-PP' 

' £i ,0.2\fn 


(.kd ) 


(-l) m+1 (±i)^-* 2 
y/nh{h + l)(2((2 + l) 


X 


^ 1+^2 

E 


„PP’ -yliMf 

L '£i /2 —m,m, 0 


kt>(kd) 


f'=Ki-A| 


(A.4) 


where the initial and final polarizations P and P’, re¬ 
spectively, may be equal or different. The ± sign has 
to be interpreted as positive if the translation occurs in 
the sense of the z axis while the negative sign applies for 
translations in the opposite direction. 

If the polarization is conserved, P = P’, the coefficient 


appearing in (A.4) reads 


<£= 2VW+i[h(h + 1 ) + £ 2 (£ 2 + 1 ) - me + 1 )] 

( A - 5) 

while for different polarizations we obtain 


cZw;m = ±4vW + 1) mkd ( P ^ P ’) • ( A - 6 ) 

In the latter case, the matrix elements of the translation 
operator vanish in the limit of vanishing wave number k 
as well as for to = 0 . 

The Gaunt coefficients (see e.g. [48] ) 


1 —m,m, 0 


(2^ + l)(2£ 2 + l)(2f + 1) 


47T 


f £1 

ll 


f h h 

e\ 

[ 0 

0 

0 J 

\ —m m 

0 J 


(A.7) 


do not require the explicit evaluation of the two Wigner 
3j symbols but can be determined efficiently by means of 
recurrence relations [55] , 

In the dipole limit, = £ 2 = 1, the matrix elements 
of the translation operator can be written as 


T-PP T-PP 
' 1 , 1 ; 0 ' 1 , 1;0 


r PP r PP 
' 1 , 1;1 ' 1 , 1;1 


t-pp't-p'p 
' 1 , 1;1 ' 1 , 1;1 



+ (kd) 3 ) ’ 
(kd) 2 + (kd) 3 ) 



2 


(A.8) 

(A.9) 

(A.10) 


These expressions make the suppression 
numbers mentioned in Sec. Ill explicit. 
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